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I. INTRODUCTION 

It is well known that topological defects may appear whenever, in the thermal history of the universe, a symmetry 
breaking phase transition occurs as for instance in grand-unified theories ffl or in some extensions of the 

standard electro- weak model [pi. Such defects represent spacetime positions where the underlying order parameter 
cannot relax, because of topological constraints, to its low energy vacuum state They are expected to interact 

mainly gravitationally with the ordinary matter so that they can induce (i) deflection and redshifting of massless 
particles, (ii) accretion of massive non-relativistic particles and (iii) emission of gravitational waves (see e.g. Q-^] for 
a review of these effects). 

One very interesting example, from a high energy physics point of view as well as from a cosmological point of 
view, corresponds to the case of cosmic strings. In this case, their phenomenological properties are determined by the 
energy density per unit length of a string, U. For instance the deflection angle is of order AttGU, G being the Newton 
constant. For defects formed at the grand-unification scale (Tgut ~ 10 16 Gev), we expect effects of a magnitude of 
GU ~ 10~ 6 . This corresponds for instance to the magnitude of the cosmic microwave background (CMB) anisotropics 
induced through the Kaiser-Stebbins effects Although topological defects may have clear observational signature 
on the CMB sky |n]|^] , observations seem to disfavor such an origin j7|-^0[ . 

Nonetheless, it does not mean that topological defects do not exist. Their detection would be of dramatic importance 
both for astronomy and particle physics Jl|,^,^ 13 since for instance estimation or bounds on their density will help 



constraining the high energy physics theories predicting their existence. Definitive predictions for string properties 
are however difficult to obtain, because in particular of the complex evolution equations that may depend on their 
microscopic structure through their equation of state. For the so-called Goto-Nambu strings (where the energy per 
unit length U and the tension T of the string are equal), it was shown that the spacetime around such straight cosmic 
string was conical pj-|l8j . Such a cosmic string formed at GUT scale would therefore induce image pairs of distant 
objects with angular separation of order 9 ~ 5.2" x (GJ7/10~ 6 ). From a pure phenomenological point of view such a 
string is expected to produce lines of double images Jl9[ . Recognizing the peculiarity of such a system, it was later 
extended to straight cosmic string with different equation of state f2(i| , pi"| and to a string with a lightlike current pulse 
PH . Moreover numerical simulations for Goto-Nambu string were also performed in the case of long strings J23] and 
cosmic loops ^4[^5). More quantitatively the prospects for a direct detection of relic string via gravitational lensing, 
and in particular the expected number of events, was first discussed by Hindmarsh Jl9| who estimated the angular 
length of string per unit area on the sky out to a redshift z to be of order 

Sloops ~ O.l^z 2 deg -1 ^longstring ~ O.lAz 2 deg -1 

where v and A are two coefficients of order unity (see also [p3[). In conclusion, it is widely believed that the observation 
of a cosmic string can be achieved through double image detections, although, in practice, it might be difficult to be 
positive about such a detection since pairs with the same angular separation appearing by pure coincidence can be 
very high (as pointed by Cowie and Hu who reported for such a cosmic string lens candidate |26|,p7|] ) . 
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The aim of this article and its companion |2£| is to make a systematic investigation of the gravitational lensing 
effects by cosmic strings. We focus our analysis on the deformation equation of a geodesic bundle in presence of cosmic 
strings (§ |l|). Standard approximations of thegravitational lens theory are also discussed in this section. After a 
description of the cosmic string dynamics in § III, we show in § [V that the deflecting potential of a cosmic string is 



we show in 

equivalent to the one by a static distribution of matter on the projection of the string worldsheet onto the observer 
past light cone. In the course in this calculation we show that the deformation field induced by cosmic strings has a 
zero convergence (without any approximation) . Examples illustrating these results are discussed as well as the validity 
of the thin lens approximation and the influence of the equation of state. In § ^ we investigate the phenomenological 
consequences of the zero convergence property on multiple image systems. In a companion paper, p8[ , we propose a 
phenomenological model of string energy distribution that gives a more quantitative account of these results. 



II. EVOLUTION OF A LIGHT BEAM 



In lens systems that are usually considered in cosmology, such as galaxies or galaxy clusters, the metric perturbations 
correspond to those of scalar perturbations. This is not the case for cosmic string effects where rclativistic motions, non 
trivial equation of state, also induce vector and tensor perturbations. We are thus forced to consider the deformation 
equations of light beams in their full generality. In the geometric optic approximation, an electromagnetic plane wave 
propagating in an arbitrary spacetime M. without interaction with matter can be described by a null geodesic [ p9[ . 
The goal of this section is to review the description of the evolution and distortion of a bundle of null geodesies and 
we start by introducing the standard elements of the gravitational lensing theory and then apply them to a perturbed 
spacetime. We then discuss the thin lens approximation and finish by some comments on its applicability. 



A. Basics of gravitational lensing 

We consider a bundle of null geodesies g propagating in a spacetime M. . Each geodesic can be described as 

g: a^(A)=i"(A)+^(A) (1) 

where (A) is the equation of a geodesic go chosen as reference and £ M is a displacement vector labeling the other 
geodesies with respect to go. Greek indices run from to 3 and A is an affine parameter along the geodesic go. With 
these notations, we can define the tangent vector to go by 

It is a null vector satisfying the geodesic equation, i.e solution of 

k^ = 0, *"V M fc" = 0, (3) 

where V M is the covariant derivative associated to the metric g^ v the signature of which will be chosen as (— , +, +, +). 

Now, we consider such a bundle converging at a point O £ M. where we assume that there is an observer with 
4-velocity u M (u M is a timelike vector, i.e. such that u^u^ — —1) and we choose the affine parameter A to vanish at O 
and to increase toward the past. We then consider at O the basis (fc M , u M , n^, ) where 2 are two spacelike vectors 
(n^nafj, = +1, a = 1, 2) such that 

<n 2/ x = 0, fc M n£ = and u M n£ = (4) 

and k^ is the null vector defined in (0). Starting from this basis at O, we construct such a basis at any point of the 
light ray worldline x^ by parallelly transporting it as 

k^V^X" = (5) 

for X = u, rii and n 2 . Note that since fc M satisfies (^) this implies that (jjj) is in fact the Fermi- Walker transport 
and thus u, n\ and n 2 remain orthonormal and satisfy (Q) for all A. Since the tangent vector k£ = k^ + d^/dX to 
each geodesic g of the bundle is a null vector, we deduce from k^k gii = g^ Jx a + £, a )k^kg — that 2fc M d£^/dA + 
k^WZ^dag^v — at first order in £. It can then be concluded that, using (||), k^ is constant along the geodesic 
and vanishes at O so that it can be decomposed as 
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a=1.2 



£o does not vanish in general, but two such decompositions (Q) with different £o parameterize the same light ray. We 
can for instance impose that £ M is spatial for the observer with 4-velocity (i.e. k^u* 1 = 0) which then fixes the 
value of £o . We also decompose the coordinates of every event of M. in the neighborhood of go as 



Xk^ + x a nZ+Tu^. (7) 



x' 

a=l,2 

The equation of evolution of £ p is obtained by writing the geodesic deviation equation |^9| 

^=K aP ^e (8) 

where R^ap is the Riemann tensor of the metric <7 M „ and where D/dA = k v \J v . Inserting the decomposition (|^) in 
(||) and using the fact that £ a = is a scalar (so that D£ a /dA = d£ a /dA with d/dA = k^d^) leads to 

I = n% (9) 

where lZ a b = R ^va(ik v k a n^n^ and a dot refers to a derivation with respect to A. Due to the linearity of the geodesic 
deviation equation (^) , £ a can be related to its initial value £ a (0) through a linear transformation T) a b as 

&(A)=Z£(A)&(0). (10) 

Since £(0) = for a bundle converging at O, with two derivatives (|Io| ) and using this equation again to eliminate 
£(0), we obtain the equation of evolution for T> a j, 

V ab = KVcb (11) 

with initial conditions 

V ab {0) = and V a b{0)=lab, (12) 

I a b being the 2x2 identity matrix. This equation has been derived in e.g. |30| |3l| . T> a i, characterizes the deformation 
field while looking in different directions. Quoting that £(0) = 9\ is the vectorial angle of observation and £(As) = As6*s 
where 9s is the vectorial angular position of the source, equation (10) can be rewritten in terms of these angles (see 
figure as 

0g = ?lMol (13) 
As 

The amplification matrix A^ = d#g/d6>j can be expressed in terms of V as 

At = BM. (14) 
As 

In the following, we decompose it in terms of convergence k and shear 7 = (p/x, 72) as 

Aab=( , \ V (15) 

V 72 1 - K + 71 J 
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FIG. 1. Description of the lensing geometry. S is the source, L the lens and I the image of S. O'S = £(As). 

B. Application to a perturbed spacetime 

We now restrict our study to a perturbed spacetime of metric 

with rj^ v being the Minkowski metric rj^ — diag(— , +, +, +). We work in harmonic gauge 

d v U" v = 0, 

so that the Ricci tensor at linear order in the perturbation h^ v reduces to 
A being the Laplacian A = did 1 , Latin indices running from 1 to 3. The Einstein equations take the simple form 



(16) 
(17) 

(18) 
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(df - A) V = 16^G(T M „ - -v»vT£) = mnGT^, 

T^ v being the stress-energy tensor of the matter perturbation. 

The solution of this equation can be obtained by means of the Green functions, of the d'Alembertian 



1 6{t' -t±\x-x' | 



47T \X — X' 

so that, using the retarded solution, we solve the Einstein equations ( |l8| ) at linear order as 

d 3 x' 



h lxv {x,t)=AG J \ x >_g\ - F , lv (x',t- \x~x'\). 
Now, we can solve (^lj) order by order: at zeroth order, lZ a b = so that it reduces trivially to 

at first order, the equation of evolution ([ll]) gives 

^i 1 6 ) (A) = A^(A) 
from which we deduce that T>^ (A) is given by 



(20) 



(21) 



(22) 



(23) 



= / " A(A S - X)U^(X)dX. (24) 



(i 



It is interesting to note that while T> a b is not symmetric in general, it is symmetric at first order in the perturbation. 
Using the expression of the Riemann tensor as 2R^ aup = h aUjiip + h^ p<crv — h vp , ap — h ap ^ p , we get that 



K$ = 7;K^ P k»k°n» a n p b - -— (T a ppVap k^n^ b ) (25) 



where k^, and are evaluated on the unperturbed geodesic (and are thus constant) and where T"^ are the 
Christoffel symbols at first order in the perturbation. Choosing 



< = K (26) 

and defining the deflecting potential $ as 



= -V*"*". (27) 



where t(X) = to — xtt (A) is the equation of the photon trajectory [to is the time of the observation to = t(X = 0)], we 
obtain that 

V^(x a ,X s ) = (j^ b A(A s -A)a afc $(^,x||(A),<(A))dA^ + jf & A(A S - A)* ab (x a , x {{ (A), i(A)) dA (28) 

where 

^ ab (S,t) = -~(T^rj ab ^) (29) 
and where d a refers to a derivative with respect to the coordinates x a as defined in (^). 

C. The thin lens approximation for static distribution of matter 

In the thin lens approximation, one assumes that the effect of the deflecting body takes place over only a small 
fraction of the light path. This approximation is usually considered in cases of scalar perturbations. The aim of this 
paragraph is to recall its derivation in the standard case to see to which extent it applies for cosmic strings. We thus 
assume that the lens is localised at A = Al with an extension small compared to Al and that this matter distribution 
is static so that T^k^k 11 = H(x±,x\\(Xi,))6(x\\(X) — arn (Al)), where E is the surface energy density. It follows that 
the deflecting potential reduces, after integration over the direction of propagation, to 

*(2 ±) x\\)=2G [ = E(xl,Z||(A L )) (30) 

where x± = (x%,X2). Since only d a b$> enters the expression of and since this quantity varies as [xn — :e||(Al))~ 3 
as soon as we are looking close to the string [i.e. when \x^ — (^|| — ^(Al))] and we can approximate the 

deflecting potential as 

$(x) = $(x ± )5(a;||-X||(A L )) (31) 

with 

$(£]_) = f S $(f ± ,a;||(A))dA 
Jo 

= 2G J E(fl,a;||(A L )) [in (a:,, (A) - x {] (A L ) + yj(x ± - x' ± ) 2 + ( X]] (A) - x\\ (Al)) 2 )] ^ d 2 x' ± . (32) 

Now, if we assume that the impact parameter is small compared to the two distances lens-object and observer-lens, 
i.e. 
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W± - S±,h\ < (ac|| (As) - ^||(Al),^||(Al)), 
we deduce that the deflecting potential integrated along the line of sight is given by 

$(£l) = -AG { ln\x ± - x / x \Jl(x' x ,xu(\ L ))d 2 x* x 



(33) 



(34) 

up to a constant which depends only on x» (Al) and x\\ (As); we forget this constant since T> a b involving only derivatives 

of $ and is thus independent of its value. The second contribution of T>^J involves the computation of the potential 
^afc and one can show from ( *29| ) that if we deal only with scalar perturbations (i.e. such that hoo = 2<p and hij = 2ip5ij) 
then = 0. For the vector and tensor perturbations, ^ a b does not vanish but in the thin lens approximation its 
contribution corresponds to a boundary term (time dependent but identical for all light rays joining the source and 
the observer) which can thus be dropped. 

Then, the amplification matrix, in the thin lens approximation, reduces to 



A a b — lab — 4G- 



. Ah — At 



As 



'-Xhd a d b / In \x± 



T,(x ± )d 2 Xj 



which rewrites, after the change of variables 8' = x^_/As, as 

As — Al 



A a b(Xs 



In 



AG- 



Ac 







f ln 


0*1-0' 



A L S(0')d 2 0' 



(35) 



(36) 



Decomposing £(0~) as AlS(0) = J (J,(s)5 (j) ' — s tring(s)) ds where s trmg represents the locus of the string on the 
plane xn = x\\ (Al), we get that (*3(|) reduces to 



A a b(Xs) = hb - d e? dg b cp(0i) with ^(0i) = AG 



As — Al 
As 



In 



0string(s) /J-(s)ds 



(37) 



where /i(s) is the projected total lineic energy density of the string (which mixes the effect of the lineic energy, the 
tension and the currents along the string if any) and ip is the effective projected potential. 

When dealing with topological defects, there are different reasons why such an approximation may not hold. 
First the strings are extended and move with relativistic speed so that (i) they are a priori not confined in a plane 
Al ~ constant and (ii) one cannot assume that the distribution of matter of the lens is static so that the time 
dependence in the line-of-sight integration in ([p]) has to be taken into account. These issues will be addressed in 



IV. B after a description of the general stress-energy tensor of strings (§ [II). 



D. Comments 



1. Gravitational potential and deflecting potential 

As a first comment, let us stress that in general the deflecting potential $ is different from the Newtonian gravita- 
tional potential. For instance, a general perturbed spacetime has the general post-Newtonian metric 

ds 2 = -(1 - 2cj})dt 2 + 2A i dx i dt + [(1 + 2i/j)6ij + 2E lj ] dx l dx 3 (38) 

where <f> and ip are the Newtonian potentials, Ai and Eij are the vector (rotational) and tensor (gravitational waves) 
perturbations satisfying 

E\ = diE ij = d t A l = 0. (39) 

ft follows that 

§ = <f> + i/) + A i i i + Eij 7 V (40) 

where j l is the direction of observation. This includes effects from the rotation of the deflecting body and of gravita- 
tional waves. Indeed, in the case of pure scalar perturbations, we recover that 

$ = 2<f). 

In the case of scalar perturbations, one can easily check that ^> a b = but topological defects also generate vector 
and tensor perturbations. In the thin lens approximation, the contribution of these two terms reduces to a boundary 
term that can be neglected but in the general case of extended object, one has to check that it is still the case for the 
vector and tensor modes. 
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2. Deflection angle 



In a general spacetime, the deflection is not straightforward to define. This is for instance the case in a perturbed 
spacetime with perturbations on all scales (see e.g. H||6| for a discussion and a generalization of this concept). If 
the matter perturbation causing the lensing is localized in space then the metric perturbations generally die away and 
the spacetime is asymptotically unperturbed. In that case, one can compute the deflection angle simply by solving 
the geodesic equation (|^) at first order in the perturbations. For that purpose, we decompose the tangent vector to 
the geodesic g as 

k^kf + Sk 1 *. (41) 

is the tangent vector of the light ray in the unperturbed Minkowski spacetime . Note that fc M is different from 
the vector fc M defined in (|J) which labels the geodesic of reference go- At first order in the perturbation, it is in fact 
possible to choose the unperturbed geodesic as reference since the displacement £ is first order in the perturbation. 

Since both rj^k^k" = and g^k^k 1 " — 0, we deduce that k^Sk^ = at first order in the perturbations. At linear 
order, the geodesic equation ([|) implies that 

^-8k a + 5T a k»k v = (42) 
clA 

with ST" = \r] af3 (hp^ v + hp v<fl - h^p), so that 



[Sk a ]° Xs = - ? r [V^]° As + \v ap I W^^dA, (43) 

•'As 

the integral being performed along the unperturbed geodesic. Forgetting the boundary term (which is a time dependent 
term but identical for all light rays joining the source and the observer), we can extract the variation of the the photon 
energy measured by an observer with velocity u M 

SE= f dt${x_L,x\\(\),t(X))d\ (44) 

and the deflection a in the plane perpendicular to the line of sight 

a = V± [ $(xl,a;||(A),t(A))dA. (45) 

The effect on the photon energy, and thus on its redshift, is nothing else but the well known Sachs- Wolfe effect j37j. 
Focusing on fl4q), since $ is evaluated along the photon geodesic t(A) = to — xii (A), we deduce that 

[d t + d Xll ] $(£]_, an (A), t(A)) =0 (46) 

along the photon path. Hence, rewriting the three dimensional Laplacian as A = 9| + Aj^ where Aj^ is the two 
dimensional Laplacian and using the Einstein equations (|l^), we deduce that 

V±.a = -8nG T (x±, x\\ (A), t(A)) dA = -8nG / T (x ± , x {l (t), t) dt (47) 

J\ s Jo 

after choosing the parameterization A = to — t and where we have introduced T = T^k^W = T^k^k" . It follows 
that the 2-divergence of the deflection depends only on the projection of T^v onto the photon trajectory in between 
the source and the observer and thus vanishes on all directions which do not intersect the string worldsheet. Note 
that such a result holds for any relativistic and/or extended lens. 

Now, in the thin lens approximation, one can relate #s and 9\ (see figure |l|) by 

61= Of - As ~ Al q q - (48) 
As 

The amplification matrix being given by A% = d#g/d0j , we obtain that 



o 



7 



A ab = I ab + As A Xh d eb d a [ " $dA. (49) 



o 



Since the lens is localised in planes close to A = Al, we have that dga w Al<9 q so that the former expression of the 
amplification matrix reduces to (^) once we use the expression ([mJ) for the integrated potential. In both approaches, 
we find that the deflection angle is given by the usual expression (see e.g. Q]) 

a = 4GAl / i ~t Z(6')d 2 6' 

J \e l -9'\ 2 

and we emphasize that defining the amplification matrix through the deflection angle implicitely assumes that we are 
in the thin lens approximation (as seen for instance on equation (Hq)). 



3. Applicability to cosmology 

It has probably not escaped a careful reader that we have restricted our calculations to perturbations around a 
Minkowski spacetime. The justification of such a choice is that the null geodesies of two conformal spacetimes are 
identical, so that lensing effects are the same. The only difference when considering an expanding Friedman-Lemaitre 
universe will be the computation of the distances, i.e. of As and Al in (36). Note also that the expansion of the 
universe affects the dynamics of the topological defects network (see e.g. [38]) but this will not be relevant on time 
scales of order of the impact parameter which are small compared to the dynamical scales of the universe. 

III. DYNAMICS OF COSMIC STRINGS 

The determination of the amplification matrix A requires the knowledge of stress-energy tensor of the cosmic 
strings. In this section, we first present the definition of this tensor and derive the equation of motion of a string. We 
then focus on the particular case of a non-rotating cosmic loop. 

A. Equations of motion of the string 



As shown by Carter p9| , there is an elegant way to describe the dynamics of a [d < n)-brane embedded in a n- 
dimensional spacetime. In this section, we recall the main steps of this formalism necessary to obtain the dynamical 
equation of evolution of the string; details can be found in 
the string worldshcct 



It requires the introduction of the induced metric on 



Jab = 9» v x,a xV ,b-> ( 50 ) 

where A,B... refers to coordinates on the worldsheet and from which one can construct the fundamental tensor ff 111 
and the orthogonal projector _L A " / as 

*r = / AB <A<B. and ^=s^-e (5i) 

The covariant derivative V defines a tangentially projected differentiation operator 

V^^V,. (52) 

The second fundamental tensor is defined by 

K^" = CV M r£, (53) 

and the condition that the worldsheet is integrable, i.e. that all its elements mesh to form a well defined d-surface, 
is expressed by 

K^f - 0, (54) 

which is a geometric identity for the worldsheet. 
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In the case of a string (d = 2), the Lagrangian density £ can be expressed as 

V7d 2 C^ (4) [*" - x»(t A )] (55) 

where / is the determinant of the metric Jab defined in (pfj|). £ is distributional and not confined to the string 
worldsheet whereas £ is locally regular but confined on the string worldshcct. The string action can then be written 
either in terms of £ or of £ as 

String = J & 2 Q^ffC = J & A x^—g£. (56) 

Now, making an infinitesimal variation Sg^ of the action ( p6| ) provides a definition of the "surface" stress energy 
tensor density T^ v (confined and regular) and of the stress energy tensor density T^ v (distributional and unconfmed) 
by 

28S stling = J & 2 C^]f^8g^ = J ^xJ—gT^b 9ilv . 

These two stress-energy tensor are related by 

= -j= J d a C-/fT'"'* (4) [xf - a^(C A )] ■ (57) 

The internal coordinate stress energy tensor, T , has been projected onto a corresponding background stress energy 
tensor, T Ml/ , as in (|l]). One can then show that the general form of the dynamical equation of the string is 

V M T^ = r (58) 

f v being the force exerted on the string by any background field such as e.g. an electromagnetic field. This dynamical 
equation fl58| ) is equivalent to the more natural equation of conservation V AI T'" y = f v with f v related to f as in (p7|). 
For a string of energy per unit length U and of tension T, the surface stress energy tensor density is of the form 

= UvPu" - JW, (59) 

where u M and are respectively a timelike (u^u^ = —1) and a spacelike (v^v^ = +1) unit vector tangent to the 
string worldsheet (i.e. u v =-L£ d" = 0) so that ?7 Miy = —u^"U v + v^v v ' . The dynamical equation governing the 
evolution of the string is given by the tangential projection of (|5^ ) which, in the free case we are considering, leads to 

VW^T = 0. (60) 

This equation of evolution (pOj ) can then be solved once we are given an equation of state, i.e. U(T). Such equations are 
pro vided once the microscopic structure of the string is described. The most well known are the Goto-Nambu strings 
po| (U = T) and the Nielsen-Olesen strings [fll] (U + T = const.) and some have been obtained for superconducting 
strings Q. 

B. Application to a non rotating cosmic string loop 

In the case of a non rotating circular loop of radius R, we work in cylindrical coordinates (t, r, 8, z) and assume 
that it is lying in the plane z = z s . Parametrising the loop worldshcct as 

*ioop = *, Hoop = r± = R(t) (cos 6, sin 6), x 3 = z loop - z s = 0, (61) 

one can show that the unit spacelike vector tangent to the string worldsheet is = 6^ and we have 

V = 7^ + 7^, e^ = R8% (62) 

with 7 = (1 — R 2 )^ 1 / 2 being the Lorenz factor associated with the radial contraction/expansion of the string. They 
satisfy _L(^ u v =-L£ 6 V = 0. From ( p^ ) and (p7|), the stress-energy tensor entering the Einstein equations ( [l9| ) is given 
by 
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T^{x,t) 



( ~f 2 U j 2 RU 

j 2 RU -f 2 R 2 U 


\ o o 










-R 2 T 
/ 



8{z - z s )5(r - R). 



(63) 



We now need to write the dynamical evolution equation (|60|) to get the evolution of the radius of the loop as a function 
of time. Using that fj^ u — —u^u v + 6^6 v and the expression (j5TJ) of the stress-energy tensor of the string, the equation 
( |60| ) can be rewritten as 



(-Uftu" + M V )V„ {Uu»u p - TO^O") = 0, 



which reduces to 



df 



(jUR) = 

^Lr = -J-1 

dt 2 j 2 RU' 



(64) 

(65) 
(66) 



This system of equations for (U, T, R) is not closed and can be solved when one specifies an equation of state U(T). 
Equation (pq) shows that the total energy jRU of the loop is a constant of motion. Note that we have not used the 
identity (d4|) which is identically satisfied in our present example. 



IV. LENSING BY A COSMIC STRING 



A. A first example 

As a first application, let us consider a static straight cosmic string lying along the axis x 2 in a plane perpendicular 
to the line of sight (direction x 3 on figure with (p — 0) so that 




«(Al0i)*(s|| -3||(Al)), (67) 

with U and T depending on x 2 — \^0 2 only- Since T — XI {x%)b{x\)&(xn — xm(Al)), the first integral of (|28|), after 
integration over x[ , reduces to 

G ( 

A ab = I ab + 2 — / U{x' 2 )dx' 2 d ab J{{x 2 - x' 2 ) 2 + x\) 

with 

J((x 2 - x' 2 f + x\) = f S ,^ As ~\ F^f dA 

Jo y/(x 2 - x! 2 ) 2 + x\ + (A - A L ) 2 

where we have chosen a parametrisation such that xu (A) = A. This latter integral can be computed and gives 



J(A) = i(As - 3A L )v/A+(A s -A L ) 2 - (\ s - ^A L \ 



+ (As A L - A 2 -\ A) In ^±^A^^1 
2 -A L + ^A + XI 

with A = (x 2 — x' 2 ) 2 + x\. It follows that the amplification matrix is given by 

A ab = lab + 2 Y s l U ( X '^ dX 'l ( 2J '( A ) S *b + {Xa ~ x' a ){x b ~ x' b )j"(A)) (68) 

where J' = dJ/dA and x\ = 0. Now, one can estimate the dominant term in the integral of expression ([38|) when we 
are looking close to the string (i.e. when x\,x 2 <C As — Al, Al). For that purpose we assume that As — Al ~ Al and 
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set A ~ As — Al — Al- We then split the integral of expression (pq) in a contribution where A 



< A 2 and another 



where A 



> A . Using the expansion of J (A) in these two regimes as 



J (A) - - (A s - A L ) 




constant + O 



A/X 2 < 1 
A/X 2 > 1 



one can show that, as long as we are looking close to the string, the dominant contribution to the amplification matrix 
is given by 



Aab 



4G- 



,As — At 



As 



d g *d e b 



In 



U{9' 2 )&6' 2 . 



(69) 



In the particular case where U is constant, this approximate give the general result of the deflection by a straight 
cosmic string and one can thus thing that this is a good approximation when U is fluctuating around a mean 

value. Note also that for such an infinite string lying in a plane perpendicular to the line of sight we recover the 
general form ( |37| ) of the deformation matrix in the thin lens approximation. 

In the following of this article, we investigate more general results concerning the deformation fiels by a cosmic 
string which do not assume that we arc in the thin lens regime. 



B. General Results 



In the general case, the source term generated by a cosmic string will be localized on the string worldsheet so that 
[see equation j59|)] 

F{x, t) = J dCF(£, t)8 (x - r(C, t)) (70) 

where (t,r(£,t)) is a parameterization of the string worldsheet; r(£, t) represents the locus of the string on each 
constant time hypersurface and T is the energy density per unit length [note that we have chosen a parametrisation 
such that t is both the coordinate time and an intrinsic coordinate of the string worldsheet which implies that we 
have a one dimensional integration on the spatial internal coordinate £ and not a two dimensional integration as in 
(|59|)]. Inserting this decomposition in we deduce that the deflecting potential ( |27| ) is given, after integration over 
space, by 

= 2G J _ d *^ f ? [r{U)A S(t'-t + \S-r(C,t')\) . (71) 
Following GjJ, the integration over t' can be performed by introducing i s tring(£, C t) solution of 

*-*stringO?,C>*) = \X -r(C, Wing C> *)) I ( 72 ) 

so that 

from which we deduce that 

$(x,t) = 2G [ - % tst ""gf ' C ' *)). *Btring (g, C, *)] dC (73) 

J \x - r(C, i s tring(a;, C, *))| ~ o t r{C, i s trin g (a;, C, t)).(x - r((, t stling (x, C, t))) 

$ on the point (x,t) is then given by the projection of the string energy on the past light cone of this point, i.e. on 
the curve {r(£, 4tring(^i Ci £))i_£string(a?, Ci t)} which is the intersection of the string worldsheet with the past light cone 
of the event {x, t) [see figure |j . 

Now, focusing on k, the deformation matrix is explicitly given by (|2^) and ( p8[ ) and, proceeding as for the deflection 
angle, one can easily sort out that 
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S ab V%= f " A(As-A)A x $(?L,a||(A),t(A))dA+ / " A(A S - A)rt afc dA 
Jo Jo 

" A(A S - A) (d 2 t - dfj ®(x±, x\\ (A), t(A)) - 87rGJ c '(fj_, a;,, (A), t(X)) 



dA 



As 



+ / A(A S - A)£ ab * afc dA 



(74) 



The first term of the first integral vanishes when evaluated on the photon geodesic so that the contribution of the 
first integral to the convergence k defined in (jig) reduces to 



k(x±, to) — AttG 



A(A S - A) 



J 7 (xj_,xi|(A),i(A))dA = 4ttG 
A3 j* 



t 



(t-t, 



emission 



)(to-t) 



^0 — t 



emission 



J r (f_ L ,X||(i),t)dt. (75) 



£f)] = 



where we recall that to = t(X = 0) is the time of reception and where iemission = t(Xs) is the time of emission. The 
contribution of the second integral of (f74|) vanishes since (i) ^ a bS ab = both for scalar perturbations [4> and ip in 
and for vector perturbations [A4 in (|38|)| and (ii) for tensor modes [jEy in (||||)] , 5 Q { ) £ afc oc ^ [(^* + ^||)(-^i 
when evaluated on the photon trajectory. We conclude that the convergence K is given by ([75|). It is then given by 
the distribution of matter evaluated on the photon trajectory, up to a geometrical factor. The lensing by a cosmic 
string is thus equivalent to the lensing by a linear distribution of matter. As a consequence, k = everywhere but 
on directions such that the observer past light cone intersects the string worldsheet; this result, valid whatever the 
equation of state, is one of the main results of this article. It holds for any relativistic lens and does not rely on the 
thin lens approximation. 

For instance if the string is lying in a plane perpendicular to the line of sight then (|75|) reduces to 



k(x±, to) = 47rG 



As — Ai 
As 



AlS(xj_,2;||(Al),<(Al)) oc Vj_.a. 



(76) 



Under this form, again we see that k = everywhere but on directions intersecting the string worldsheet. 

Note that the drivation of ( |7q ) relies strongly on the fact that we took the trace of the deformation matrix in < \74 
Therefore, similar expression to (|T5j) cannot be obtained for the other components of the amplification matrix and 
one has to rely on ([73h. 



It has also to be noted that the expression of the amplification matrix in (49) (which relies on the thin lens 
approximation) together with the general result ( |75| ) implies that the phenomenological description of A% in ( |37j ) is 
very general. It holds for any string dynamics provided the extension of the string is small enough for the thin lens 
approximation to hold. Different aspects of these results are illustrated in the next paragraphs. A more elaborate 
phenomenological investigation based on Eq. d37]) is proposed in a companion paper pq| . 




FIG. 2. The intersection of the loop worldsheet and the past light cone of the event (x,t). The dash circles represent the 
loop in different constant time hypersurfaces and the dash-dot lines the loop worldsheet. 



C. Lensing by a non— rotating cosmic string loop perpendicular to the line of sight 

We now consider as an application the case of a non-rotating circular loop oscillating in a plane perpendicular to 
the line of sight. Its dynamics is described by the set of equations (|65]-p6[) and from its parameterization ( |6l| ) we 
deduce that 
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F(x,t) = R(t)j(t)U(t)S(x 3 )S(x x - f x )d0. (77) 



The deflecting potential <Y21\) integrated along the line of sight is then given by 



J[$] = / S $(x(A),i(A))dA 



+00 



6 (V + X 3 - r + yf(x± -X'^f + X 2 ^j 



= 2G dx 3 d6 d 2 f x / dt'-yRU^ = ^<5 (f ± - fj_(t', (9)) (78) 

where we have parameterized the geodesic as t = tq — x 3 with tq = to ~~ x h ■ Now, defining the new variable z as 



z = x 3 + J(xx - + x\ (79) 



which satisfies 

dz dx 3 



\J(x ± - x' ± ) 2 + xj 



(80) 



The integrated potential ( |78[ ) reduces to 

f27T J _ r + oo 



/[*] = 2G / dO d 2 f ± / — / dt'-fRU5{t' + z- t )S (x' ± - r±(t', 6)) (81) 
Jo Jr. 2 J a z J~oo 

where the limits of integration are 

A = -x L + J(xx - x' ± f + x\~ \ (f± ~ ^ (82) 



2 XL 



B = xs —XL + y (S± — £\) 2 + (x s - x L ) 2 ~ 2(x s - x L ). (83) 

The approximate values of A and B are obtained at lowest order when we consider zones close to the string in 
comparison with the distance string-observer and string-source. 

After integration over t' and using the loop equation of motion (|65|-|6rj|) which state that r yRU is constant, we get 
that 

I[$] = 2G~/RUJ(x x ,t ), (84) 
where J(x±, to) is a dimensionless geometrical integral given by 

r /.2(x s -:cl) J 

J@±,to)= d9 d 2 af x / (gi _ g ,Y — S(x" ± -r ± {-z-T ,e)). (85) 
Jo Jr 2 Jh-^ 1 — — z 

This integral can be rewritten as 

J(Sx,t )= — dO d 2 x' ± 5(x' ± -rx(-z-T ,e)). (86) 



z 

We deduce that, at z constant, the integral over x' ± reduces to the computation of the angle {3 (see figure ||) of string 
within the disk of radius y/2xjjz which can be computed as followed: 
Setting u = ^/2x\~z and h = OH, we deduce from 

u = u 1 +u 2 , u\ + h 2 =x\, u 2 2 + h 2 = R 2 (87) 

where u\ = CH and u 2 = HD on figure || and R now stands for R(t = to), that 



13 



x 2 , - R 2 + u 2 

Hi = U 2 

2u 



x\ + R 2 + u 2 
2v\ 



and thus that 



_ xi+R 2 - u 
cosp = 



2\x±\R 



(89) 



The integral over 9 and x' ± obviously vanishes if u < \\xj_ \ — R\, reduces to 2n when u > \x±\ + R and gives 2j3(u) 
otherwise. Then, after splitting the integral over u in ( p6[ ) in three pieces ([0, — R\], [\\x_\_\ ~ R\ , |xj_| + R] and 
[\x±\ + R, 2 v /x L (x s - x L )]) we get 



J(f±,t )= / 2f3(u)— + 2n — , 

J||xx|-R| u J\x ± \ + R u 

with /3 given by ( |S9l ) . We can compute this integral in the two following regimes 
1. If |xj_| < R{tq), can be rewritten as 



J = 2\x±\R 
which can be computed to give 



Arccosw 



! |f_L| 2 + i? 2 -2|xj_|i?U 



dv + 2tt In 



2^x L (x s - z L ) 



J = Ci 



(90) 



(91) 



(92) 



where Ci is a constant depending on xi,, xs and R(tq). Then, there is no deflection of a light ray passing inside 
a large loop, as first pointed out in pa] and as expected from the Gauss theorem. 



2. If lafjj > R{t~o), @) now gives after integration 



J = d - 2tt In 



Fx I 
R 



(93) 



and we conclude that a small loop perpendicular to the line of sight acts as a point mass M = 2ir"fRU whatever 
its equation of state. We checked that this is also valid for a tilted circular loop and it is natural to expect that 
the fact that a loop acts as a point mass at a distance larger than its caracteristic size is valid whatever the 
geometry of the loop. 

One can also check from (|92]-p3|) that Aj_J(xj_,t ) = if xj_ ^ R{t ). Since k oc Aj_J, we recover the result 
from that the convergence vanishes if x± =^ r± . 





o>V\ 






\ ;i 




1 R 








/ pi,-' 

/*^XP2 , 
I?' K 



FIG. 3. Integration over the disk of radius u = \/2xlz around the point C = x\ 
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D. Tilted static straight cosmic string 



To finish, let us consider a tilted static straight cosmic string aligned along an axis making an angle ip with the 
direction x 2 [see figure |J and for which, from ( p7| ) , 



T = \ [U(£)-T(i) sin 2 <p] d£S(x - r(£)) 



where f(£) is a parameterization of the string. If we choose I such that 

f(i) : ^1=0, r 2 =£cosip, r 3 = X|| (A L ) — £ sin ip, 
the deflecting potential ( |73| ) is given by 

df 



$(x) = 2G 



%£)\ 



[U-T sin 2 ip] {£). 



When U and T are constant, it can be integrated to get 

$(x) = 2G [f7 - T sin 2 99] (Coo - In [x 2 - (x 2 cos ip - x 3 sin ^) 2 ] ) 



(94) 



(95) 



(96) 



(97) 



where Coo is an infinite constant and where we have introduced X3 = xu — xl- The infinite constant Coo can then 
be forgotten because only the derivatives of <E> enters the computation of the deflection angle and of the amplification 
matrix which are the observable quantities. The deflection angle is then given by 



§(xi,X2,Xz)&Xz. 



(98) 



After integration over X3, we get 

4G [U-T sin 2 ip] 



Oil 



a 2 



arctan 



cos ip 

4G [U-T sin 2 ip] tan ip 



^x 2 sin 4- cos s 



X\ 



-1 X 3 =XS— X L 



hi 



\J x\ + (X2 sin ~tp + X3 cos ip) 2 



x 3 =x s -Xi J 



(99) 



In the limit where (xs — xl) |cos<^| and xl |cosy| are much larger than |xi| and \x 2 sin^|, we get 

ai ~ 47rG [U cos ip - 
a 2 — AGU tan ip In ■ 



ai ~ 47rG [U cos + (U — T) sin ip tan y?] 
xs - x L 



(100) 



XL 

This has to be compared with the standard result for a Goto-Nambu string for which a — AttGU cos ip ||,[lf|[l9|,^0| . 
Now, as pointed out by Peter plj in the case of a string perpendicular to the line of sight, there are two origins to the 
deflection: the deficit angle (term proportional to U +T) and a contribution from the curvature (proportional to U—T). 
One can understand such a result by decomposing the stress-energy tensor ( |67| ) as 2J r A1 „ = 2U x diag(0, 1, 0, 1) + (U — 
T) x diag(l, —1, 1, —1), i.e. as the superposition of a Goto-Nambu string and a linear distribution of non-relativistic 
matter of density p=U — T per unit length. Then it is straightforward to see that the bending angle of the second 
contribution depends only on the projected mass per unit length and thus becomes larger by a factor 1/ cos ip as found 
in ( |100| ). A consequence of this result is that, for general cosmic strings not perpendicular to the line of sight, one 
expects to have larger deflection than for a Goto-Nambu string. 
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FIG. 4. Parameterization of a static straight tilted cosmic string. 



In this case, the accuracy of the thin lens approximation can be investigated. For that purpose, we remind that the 
shear is given by 



7i 

72 



= 2G 



(ar L + x 3 )(x s - x h - x 3 ) (\{dl- df] 



•i'S 



J 2 

did 2 



[U(£) - T(£) sin 2 <p] dl 



\Jx\ + (x2 — I cos (p) 2 + {x 3 — I sin ip) 2 



(101) 



Due to the derivatives with respect to x± and X2, it is easy to see that the integral over £ is peaked around x 3 
£suiip ~ xi tan^J. Thus, on a field of width xi — xs02, the variation of geometric factor is bounded by 



<5|(x L - X 3 )(x S -£l - x 3 ) 



1 + 2^ 



xl/xs (1 - x L /x s ) 



| tan v| ^2 < 



x-l/xs (1 - xl/xs) 



|tany!| 6> 2 , 



from which we deduce that since 02 <C 1, the thin lens approximation is still very good for tilted string with a tilt not 
larger than ip = n/A say (see figure wfor a numerical estimation of the relative error). 




0.2 0.3 



FIG. 5. Relative error on the geometric factor with respect to the thin lens result for a string with a tilt ip — n/4 on a field 
82 — 10". The relative error scales like (82/W") x tan (p. 



E. Discussion 



In this section, we have shown that the deflecting potential of any extended lens with relativistic motion is obtained 
by considering the projected energy density on the photon past light cone. This implies, in the case of cosmic strings, 
that the convergence n vanishes everywhere but on the string projection. 

We then studied the case of a loop oscillating in a plane perpendicular to the line of sight and show that the 
equation of motion of the loop can be used to integrate the deflecting potential. We found that a photon propagating 
inside such a circular loop was not deflected and those propagating outside were deflected as if the loop was a point 
mass object. This generalizes the result by de Laix and Vachaspati j24| to strings with any equation of state and 
shows how the equation of motion of the loop enables to factorize the integrated deflecting potential. This lets us 



16 



conjecture that this result will be valid whatever the geometry of the loop, the geometric factor J being different 
for each loop geometry. We finished by discussing the case of static tilted cosmic strings to emphasize that, for 
non-Goto-Nambu strings, there can be larger deflections and we also discussed on this example the validity of the 
thin lens approximation for strings. 



V. PHENOMENOLOGY OF A DEFORMATION FIELD WITH « = 



As seen in the previous section, we expect the deformation field of a cosmic string to be such that k — 0. Indeed, 
in ( frq ) we only showed that k vanishes in all directions such that the light ray arriving with this direction does not 
intersect the string worldsheet. In this section, we are mainly interested in objects which do not overlap the string 
such that their deformation is the one with a k = field. Another restriction of this study is that we assume that 
the caracteristic size of the object is smaller than the caracteristic size of the variaiton of the shear 7; thus we donnot 
consider galaxies lying in the immediate neighborhood of the string (see the companion article |2§| ] for an illustration). 
With these restrictions, we can consider that the deformation field has a zero convergence and the goal of this section 
is to study the main phenomenological properties of such a field. The two kinds of questions we would like to answer 
are: 

(Ql) Given a source of surface S s and ellipticity e s , what can we say about the morphologies of all its possible images? 

(Q2) Given two objects, how can we know that they are the images of the same source? This question reduces to 
study the allowed morphologies of the sources that have the same images. 

We start by setting the general framework and then study respectively (Ql) and (Q2). This study is a first step 



toward the discrimination between pairs of lensed sources by a cosmic string and fake lenses [ 
made within the assumption that the shear variations over observed background images is sma 
a severe limitation for such an approach if the string energy density is small. 



2J,|27j. This study is 
1. It may actually be 



A. Describing the morphology of a cosmic object 

To any object of elliptic shape such as a galaxy or a cluster, we can associate a positive definite symmetric matrix 
M a b describing its shape as 

X l MX < 1 (102) 
where X 1 = [x\, £2)- This matrix can be diagonalised as 

M = pW A + a ° W) (103) 



where P is a rotation matrix defined by 



cosy sm6 
w 1 — sin cos ' 



and a subscript t denotes the transposition. A_ < A + are the two positive eigenvalues of M and 9 is an angle describing 
the orientation of its principal axis with respect to the basis n£. Thus any object can be characterized by the set (8, 
A_, A + ) from which we can define the surface S and ellipticity e of the object respectively as 

S{M) = det(M) = A+A_, (105) 



and we also define e as 



^l-e 2 = 4^1. (107) 



[tr(M)] 2 ' 

These definitions can indeed be inverted to get the two eigenvalues in terms of e and S as 
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4=sQ±£f. (108) 

Following Mellier , the cllipticity must be bounded by 

e > 0.5 <=> e < 0.71. (109) 

In the following, we will not be interested in the orientation of the object and we then define the shape as being 
the set (S, e). The shape matrix M l of any image can be related to the shape matrix M° of its associated source as 
(see e.g. Q) 

M 1 = A~ 1 M S A~ 1 (110) 

[this is valid only if we consider that the carateristic size of the source is smaller than the characteristic size associated 
with the variation of the shear]. Decomposing M s as in ( |103| ) and introducing M l = PM l P l which represents the 
same shape as M 1 after a rotation of —8, we obtain that 

M 1 = (PA^P*) ( A + A ° ^ (PA^P*). (Ill) 

Thus, M 1 is the image of the source M s = f ^ ^ ^ by the transformation, 

-i D i _ 1 ( l + 7i -72 



= iM" 1 *" = - = 1T J' '£ with 7 = P(-20)t*. (112) 

1 - 7 2 V ~72 1 - 7i / 

As long as we are interested only on the shape (i.e. surface and ellipticity) of the sources and/or images, we always 
can choose one of them to be diagonal. 

B. Morphology of the images of a given source 

^From the previous analysis, we can conclude that, if we are not interested in the relative orientation of the source 
and of the image, we can just restrict the problem by considering the source and the transformation matrix to be 
given by 

"-til)- ^t^O-Jv- 7 ;)' <™> 



Setting 



7 = 7(^1 with 7 >0 and ae[0,27r], (114) 



we can easily show that the shape (S l , e 1 ) of the image is related to the one of the source (S s , e s ) by 

(1-7 ) 

fI = T, T-^T ^ eS > ( 116 ) 

(1 +7 2 + 27e s cosa) 2 



where e is defined in (107). In the case of a circular source (e s = 0) we deduce that, since (5 I /5 s ,e I ) depends only 
on 7, (i) two images of same surface have same ellipticity and that (ii) all the images lie on a curve in the plane 
(S I /S s ,e I ). In figure ||, we depict the variation of S 1 and e 1 in function of 7 and the ensemble of the images of a 
circular source. 
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FIG. 6. The variation of S 1 /S s (top panel), e 1 (middle panel) in function of 7 and the set of the images (lower panel) for a 



circular source (e 
7 > 1 (solid line). 



0) of surface S . The two branches represent respectively deformations such that 7 < 1 (dashed line) and 



In the general case (e s ^ 0) we can determine all the morphologies of the images of a given source in the plane 
(5 r /5 s , e 1 ). In figures m, we depict these sets respectively for e s = 0.25 and e s = 0.5. We note that all the curves have 
the same asymptot (5 /5 s ) = 00, e 1 = 1 which is reached when 7=1, i.e. on the critical line; on these points, the 
amplification /i = [det(y^)]" 1 is infinite. Whatever e s , there exist always two circular images (i.e. such that e 1 = 1) 
obtained for 



7± 



it 



a — n [2tt] 



(the condition on a is obtained from the requirement that 7 > 0; there are also two solutions for a 
lead to negative values of 7). Now, since from (115) S^_/S\_ = (1 — 7^) 2 /(l — 7?) 



(117) 

[2vr] but they 
the measure of the area of two 



such circular images enables us (i) to determine the ellipticity e of their common source and (ii) the shears 7± of the 



two deformations. Indeed the bound on the ellipticity (10E) has to be fulfilled by e s and this can be a test to reject 
fake pairs of images. 

In the more general case of a pair of non circular images, one cannot reconstruct the ellipticity of their source but 
we can still conclude from the ratio of their surfaces if they correspond to transformations, 71 and 72, that are both 
subcritical (71 < 1 and 72 < 1) or where one is critical and the other subcritical (71 < 1 and 72 > 1). 




FIG. 7. Set of all morpholohies ((S^e 1 ) of the image of a source of ellipticity e s = 0.25 [left] and e s = 0.5 [right]. As 
explained in the text, we see (upper right plot) that there always exist two circ ular images, one for a subcritical transformation 
(7 = 7- < 1) and one for a critical transformation (7 = 7+ > 1) [see equation(117)]. The triangles represent the images by the 
transformation with 7 = 0. 
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For small deformations (7 <C 1), we have that 

S 1 /^ ~ 1 + 2 7 2 + C( 7 4 ) (118) 
eVe s ~ 1 - 4e s cosa 7 + 4 (3(e s ) 2 cos 2 a - l) 7 2 + C( 7 3 ) (119) 

so that the images almost lie on a parabola centered on (S^e 1 ) = (S , e ). In such a weak field, two images of the 
same object will have almost the same surface but can have very different ellipticities. 



0.8 r 




FIG. 8. The images of the different sources of respective shape (1,0), (1.5,0.1), (2,0.2), (2.5,0.3), (3,0.4) for a small deformation 
< 7 < 0.3. The circles represent the image by the transformation 7 = 0. 



C. Morphology of possible sources of a given image 

We now ask the reverse question: given an image (S l , e 1 ), from which set of sources can it be the image? Following 
the description and notations of the two previous sections, we now set 

-■-(•O' "-I 1 ;*!?*) <i2o » 

from which we deduce that, since M s — AM 1 A, the shape of the source is related to the one of its images by 

5 S = (1- 7 2 ) 2 5 I (121) 

e s = i 1 -! 2 ? 1 (122] 

(I + 7 2 -2 7 e I cos a) 2 1 ' 

As a first exerc ise, we depict on figure ^| the sources (S* s , e s ) of a circular image (e 1 = 0). A priories on the ellipticity 
of the sources (|109| ) and on the strength of the deformation 7 may enable us to extract from such a plot informations 
about the source of a circular image. 

In the general case where e 1 ^ 0, we can reconstruct all the morphologies of the source that can lead to the observed 
image. As an example we depict such sets in the two cases where e 1 = 0.5 and e 1 = 0.25 respectively on figures [lfj| 
and 11. We note that all the curves pass through the point (S s , e s ) = (0, 1) which is reaches when 7=1, i.e. on the 
critical line. 

Again, one can sort out that any object can be the image of two circular sources obtained by the two transformations 

7± = gI 7=f a = 0[2^] (123) 
1 T ve 1 

and the measure of (S 1 , e 1 ) enables to determine j± and the two surfaces S±. 
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FIG. 9. Sources of a circular image. We show the variation of S s / S 1 (upper panel) and of e s (middle panel) with respect to 
the strength of the deformation 7. The two branches correspond respectively to transformations such that 7 < 1 (dash line) 
and 7 > 1 (solid line). 
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FIG. 10. Variation of the ellipticity of the source of an image of ellipticity e 1 = 0.25. The circles represent the deformation 
with 7 = 0. 
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FIG. 11. Variation of the ellipticity of the source of an image of ellipticity e 1 = 0.25. The circles represent the deformation 
with 7 = 0. 



If one measures the shape of two images (S\, e\) and (S^, e\) one c an re construct, as in figures [T(] and |Tl|, the set 
of morphologies of their possible sources. Given bounds on e s , as in (10E), and on 7 one can figure out graphically 
if these two observations are likely to be images of the same object. Indeed for very weak deformations (7 <C 1) one 
gets that two images of the same objet must have almost the same surface but can have very different ellipticity (see 
figure O where we have plotted the source shapes of objects of different shape for small deformation). 



0< y<0.3 




FIG. 12. The source shapes of different images of respective shape (1,0), (1.5,0.1), (2,0.2), (2.5,0.3), (3,0.4) for a small 
deformation < 7 < 0.3. The circles represent the deformation 7 = 0. 



VI. CONCLUSION 



We have considered the general lensing properties of objects such as cosmic strings where relativistic motions and 
non-trivial equation of state induce metric perturbations of all sorts. We demonstrated that the deformation field of 
a string system on the image plane is the same as the one of a static linear distribution of matter projected on the 
photon trajectory. A consequence of this result is that the deformation field has a vanishing convergence (k = 0) 
everywhere but on the projection of the intersection of the observer past light cone and the string worldsheet. We 
explicitly illustrate this result with the case of a circular cosmic string loop in a plane perpendicular to the line-of-sight 
for which we generalize the results found in H] to string with any equation of state. We also showed that (i) the fact 
that the deformation field outside the loop is equivalent to the one obtained by a massive point and (ii) that a light 
ray passing inside the loop is not deflected are due to the energy conservation of the string. 

We also paid attention to the validity of thin lens approximation for this unusual lens system. This approximation 
is discussed in detail through the case of a static tilted straight cosmic string. It lead us to point out that for string 
with a general equation of state, the deflection may be more important than for a Goto-Nambu string, this being 
understood by the fact that the stress-energy tensor of a general cosmic string can always be decomposed as the 
superposition of a Goto-Nambu string and a lineic distribution of non-relativistic matter. The deflection is then 
due to the combined effect of the deficit angle of the Goto-Nambu string and of the curvature induced by the lineic 
distribution of matter. 

We also studied general phenomenological consequences of deformation fields with zero convergence on multiple 
image systems, the main goal being to be able to assess if two images are likely to form an image pair of the same 
source. For that purpose, we derived all the image shapes of a given source as well as all the source shapes of a given 
image. These results may serve as a groundwork for the elaboration of string detection strategies. 

We are aware that this latter part is limited by the fact that we only took advantage of the zero-convergence 
property. In practical more intricate distortion properties of the images are likely to be useful. This is one of the 
aim of the companion paper [p8| where more quantitative phenomenological properties are presented that take into 
account local string energy fluctuations. 
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